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1 Introduction 


Ordinary incompressible Newtonian flnids are described by the Navier-Stokes 
equations. These equations have been used by engineers and physicists with 
a great deal of success and the range of their validity and applicability is well 
established. Together with other fundamental systems like the Schrodinger 
and Maxwell equations, these equations are among the most important equa¬ 
tions of mathematical physics. Nevertheless, their mathematical theory is 
incomplete and requires cut-offs. The present state of knowledge is such 
that different approximations seem to be useful for different purposes. The 
mathematical questions of existence and regularity for incompressible fluid 
equations have been discussed in many books and review articles (for instance 


[0, m, H, 0,0,0, 


^1 and many more). In this work 


I describe some results reflecting research concerning diffusive-Lagrangian 
aspects of the Navier-Stokes equations [0, 0. There are two distinct 
classes of approximations of the Navier-Stokes equations that we consider. 
In one class the energy dissipation is treated exactly, but the vorticity equa¬ 
tion is not exact. This class contains the Galerkin approximations 0, 0 
and mollihed equations |^, 0] (see (0) below). The other class treats 


the vorticity equation exactly but the energy dissipation is approximated. 
This is the class of vortex methods (see (0) below) and their general¬ 
izations. This class is related by a change of variables to a class of hltered 
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; the models [K 


are a 


approximat ions (0) of the formulation 
subclass of these. The Navier-Stokes equations and their various approxi¬ 
mations can be described in terms of near identity transformations. These 
are diffusive particle path transformations of physical space that start from 
the identity. The active velocity is obtained from the diffusive path transfor¬ 
mation and a virtual velocity using the Weber formula. The active vorticity 
is computed from the diffusive path transformation and a virtual vorticity 
using a Cauchy formula. The path transformation and the virtual helds are 
computed in Eulerian coordinates (“laboratory frame”). In the absence of 
kinematic viscosity, both the virtual velocity and the virtual vorticity are 
passively transported (“frozen in”) by the flow. In the presence of viscosity, 
these helds obey diffusion equations with coefficients that are proportional 
to the kinematic viscosity and are derived from the diffusive transformations. 
The diffusive path transformations are used for short time intervals, as long 
as the transformations do not stray too much from the identity. The du¬ 
ration of these intervals is determined by the requirement of invertibility of 
the gradient map. If and when the viscosity-induced change in the Jaco¬ 
bian reaches a pre-assigned level, one stops, and one restarts the calculation 
from the identity transformation, using as initial virtual held the previously 
computed active held. 


2 Energy Dissipation 

An incompressible huid of constant density and temperature can be described 
in terms of the huid velocity m(x, t) and pressure p(x, t), functions of Eulerian 
(laboratory) coordinates x, representing position and t, representing time. 
The Navier-Stokes equations are an expression of Newton’s second law, and 
in the absence of external sources of energy they are 

dtu + u ■ Vu — uAu + Vp = 0, (1) 

coupled with the constraint of incompressibility 

V-M = 0. (2) 

In the initial value problem, a velocity 

u{x, 0) = Uq (3) 
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is given at t = 0. The coefficient z/ > 0 is the kinematic viscosity. In an 
idealized situation the velocity is dehned for all x G and vanishes at 
inhnity. The vorticity equation is obtained by taking the curl of the Navier- 
Stokes equations: 


dtU! + u ■ Vu! — uAoj = ui ■ Vu 


(4) 


where 


a; = V X M 


(5) 


is the vorticity. If u is divergence-free one may invert this relation; Dehning 
a stream-vector ip that satishes —Aip = uj using the Newtonian potential and 
then taking its curl, one obtains the familiar Biot-Savart law 

u 

If the initial velocity vanishes then u{x,t) = 0, p{x,t) = 0 solve the 
equations. Moreover, if uq is close to Mq = 0 then the solution u{x,t), p{x,t) 
exists for all time, is smooth and converges to 0. The open question in this 
situation concerns the behavior of the solution for large initial data. The 
Navier-Stokes equations are a parabolic regularization of the Euler equations 
(obtained by setting u = 0). Although the viscous term is important, for the 
study of large data one needs to consider properties of the Euler equations. 
The Navier-Stokes equations conserve momentum (integral of velocity in the 
present setting). The total kinetic energy 


[X = - 


An 


'R3 


k - y\' 


{x-y)x uj{y)dy. 


( 6 ) 
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u{x^ t)\^dx = K(t) 


(7) 


is dissipated by viscosity 

- I \u{x,t)\‘^dx + I I iVui^x, s)\‘^dxds < - j \u{x,to)\‘^dx. (8) 

2 Jb.3 JrS 2 

The dissipation of kinetic energy is the strongest source of quantitative infor¬ 
mation about the Navier-Stokes equations that is presently known for all so¬ 
lutions. This dissipation is used to construct Leray weak solutions with hnite 
kinetic energy that exist for all time, u G L°°{dt] L'^{dx)), Vu G L‘^{dt ® dx) 
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This class of solutions is very wide. The solutions have partial reg¬ 
ularity PI but are not known to be smooth. The uniqueness of the Leray 
weak solutions is not known. Neveretheless, immediately after inception, at 
positive times arbitrarily close to the initial time, the Leray solutions have 
square integrable gradients and the solutions become smooth for an interval 
of time. The solution is then uniquely determined and remains smooth for 
an interval of time whose duration is bounded below by a non-zero constant. 
The issue is whether the smooth behavior continues for all time. The simplest 
self-similar blow up ansatz of Leray has been ruled out |^, |^. The most 


important task is to obtain good a priori bounds for smooth solutions of the 
Navier-Stokes equations. If one has good bounds then the smoothness and 
uniqueness of the solution can be shown to persist |^, |Q, [^, 


m 


In situations in which such bounds are not available, the study of solutions of 
the Navier-Stokes equations needs to be pursued by considering long-lived ap¬ 
proximate solutions. The advantage of dealing with approximations, besides 
practicality, is conceptual simplicity: one may formulate sufficient conditions 
for global regularity quantitatively, in terms of the approximate solutions. If 
one devises approximations and obtains uniform bounds for them, then, by 
removing the approximation, one obtains rigorous bounds for weak solutions 
of the Navier-Stokes equations that are valid for all time. For instance one 
can prove: 


Theorem 1. Let uq be a function in that satifies the divergence- 

free condition V ■ Mq = 0 m the sense of distributions. Let T > 0 be ar¬ 
bitrary. There exists a Leray weak solution {u{x,t), {p{x,t)) of the Navier- 
Stokes equations that is defined for t G [0,T], satisfies 


1 
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'R3 


\u{x, t)\‘^dx + h> 



to 


\'Vu{x, s)\‘^dxds < 


1 

2 



u{x, to)\‘^dx 


for all t > to and t^ E I <Z [0,T], where I is a set of full measure |J| = T. 
The initial time belongs to it, 0 G /, m other words tQ = Q is allowed. In 
addition, the solution satisfies 

u{-, s)||i,oocis < iFoo 

with Koo a length scale determined mainly by the initial kinetic energy and 
viscosity: 

Hoo = C J \uo{x)\'^dx 
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If the initial vorticity ojq = V x uq is in \ujQ\dx < oo then it remains 

bounded in If and moreover 



uj{x, t)\dx + 


1 



u{x, t)\‘^dx < 



uj{x, to)\dx + 


1 



u{x, tQ)\‘^dx 


for t > to, to G I. If the initial data is in and the initial Reynolds 
number is small then the solution is infinitely differentiable for positive time 
and converges to 0. More precisely, if 



v 


1 



\uo{x)\^dx 



I VMo(a :)\^dx 


1 3 

< 27 r 23 4 


then the solution exists for all t > 0, belongs to C“(R^) and converges to 0. 


This theorem combines the bound in |P] that was proved using a version 
of the retarded mollihcation approximation procedure of [Q with the result 
of @] that was proved in the space-periodic case using Galerkin approxi¬ 
mations. The last statement about smooth solutions is proved by studying 
the evolution of the product of energy and enstrophy. The specific con¬ 
stant (about .495) comes from the fact that, for divergence-free, zero-mean 
functions ||m|| < Vm||l 2 -^|| Ari||i 2 (see (^^) below. If the initial 

Reynolds number is small then it stays small and its rate of dissipation is 
a well-known quantity that controls global existence (see (p3[) below). Both 
the mollihcation approximation and the Galerkin truncation approximation 
procedure respect the energy dissipation inequality (|[) exactly, but they in¬ 
troduce errors in the vorticity equation. The (not retarded) mollihcation 
equation is described below. One dehnes a mollihed u by 


N<s = 




/R3 


x-y 


u{y)dy = Js{-iV)u. 


(9) 


Here 5 > 0 and the positive kernel J is normalized J{x)dx = 1, smooth 
and decays sufficiently fast at inhnity. Two canonical examples of such J 
are the Poisson kernel J{x) = 7r“^(l -f- |xp)“^ and the Gaussian J{x) = 
The Fourier transforms of J, J(^) = and, respectively 
•^(0 = exp( — ^-j, are non-negative, vanish at the origin, decay rapidly and 
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are bounded above by 1. Because of the fact that at the Fourier transform 
level one has 


= mm), ( 10 ) 

the operator of convolution with Js, [m]^ = Js{—i'V)u is a classical smoothing 
approximation of the identity. The mollihed equation is 


dtu + [u] ■ Vm — vAu + Vp = 0 (11) 

together with V ■ u = 0. Here [u] = [u]^ is computed by applying the 
molliher at each instance of time. This nonlinear partial differential equation 
has global solutions for arbitrary divergence-free initial data Uq G L^. The 
solutions are smooth on (0, T] x and, moreover the energy inequality 
is valid for any Iq G [0,T], t > to. The vorticity of the mollihed equation 
does not obey exactly (^). By contrast, classical vortex methods Q respect 
the structure of the vorticity equation (|) but do not obey exactly the energy 
dissipation inequality (|]). In this paper we call vortex methods the equations 

dtuj -\-[u]- uj — uAuj = a; ■ (12) 

with u calculated from oj using the Biot-Savart law (H), and [u] = [m]^ com¬ 
puted from u using the molliher (^ . Both equation and solutions depend on 
6 but we will keep notation light by dropping the reference to this depen¬ 
dence: UJ = uji, [u] = These vortex methods may also be described by 

using an auxilliary variable w. One considers the equation 


dtw + [u] ■ Vtc — uAw + (V[m])* tn = 0. (13) 


[M* means the transposed matrix). A direct calculation verihes that the curl 
of tc , V X w obeys the equation ([I^) , as does uj . This calculation uses only 
the fact that [u] is divergence-free. The system formed by the equation (PIB]) 
coupled with 


[u] = Js{-tV)P{w) 


is equivalent to (1^ 


Here P, 

P,, = 6ji - d,A-^di 


(14) 


(15) 
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is the Leray-Hodge projector on divergence-free vectors. The initial w is 
required to satisfy Ptco = Uq. At hxed positive 5 the solution is smooth and 
global. If 1 / = 0 then these systems have a Kelvin circulation theorem: the 
integral §^w ■ dx is conserved along closed paths 7 that are transported by 
the flow of [m]. (In contrast, the mollihed equations do not have a Kelvin 
circulation theorem). The energy dissipation principle for the vortex method 


IS 


- / u{x^t) ■ [u]{x^t)dx + V / / tr {(VM)(a;, s)(V[M](a;, s))*} < 

2 JRS Jto JRS 


2 7r3 


u{x, to) ■ [u{x, to)]dx 


(16) 


This is obtained by taking the scalar product of (|^ with [ijj] where u = 
V X Ip. One can obtain the energy dissipation principle also by taking the 
scalar product of (|T^ ) with [m]. One uses the fact that Js^—iV) is a scalar 
operator (multiple of the identity as a matrix, i.e. acts separately on each 
component of a vector) that commutes with differentiation. Then the can¬ 
cellation of the nonlinearity follows from the divergence free condition. The 
energy dissipation principle gives strong control on the mollihed quantities 
(or weak control on the unmollihed ones): 


u{x,t) ■ [u]{x,t)dx = - {J{SO) N(^,t) 

^ ./rS ^ JR3 


-1 


de (17) 


and 


/ tr {(VM)(a;, s)(V[M](a:, s))*} dx = / 

'r 3 dR3 


(18) 


Because J ^ is a positive function that grows exponentially at inhnity, the 
inequality implies real analytic control on [m]: 


IL 


{Js{-iV)) 2 ([M])(a;,t) 


dx+ 



R3 


V(J,(-zV))-nM)(a:,s) 


dxds < 7 / \uo{x)\‘^dx (19) 

2 Jus 


One needs to bear in mind, however, that this is a weaker bound than the 
bound provided by the energy dissipation (||) for the mollihed equation (phj), 
where {Js{—iV))~^ [w] is bounded in 
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3 Uniform bounds 


The energy dissipation principle (H) holds exactly for the mollified equation 
(1TT|) and has a counterpart for the vortex method (|T^ in ([l6| , |^. These 
are uniform inequalities, in the sense that the coefficients are ^-independent 
and the right hand sides are bounded uniformly for all 5 > 0. Most uniform 
bounds are inherited by the solution of the Navier-Stokes equations by pas¬ 
sage to limit. Some uniform bounds for the equation (|TI]) can be summarized 
as follows: 


Theorem 2. Letuo be a square-integrahle, divergence-free function. Let 6 > 
0. Then there exits a unique solution {u,p) of defined for allt > 0. The 
solution is real analytic for positive times. The limit lim^^o= Uo{x) 
holds in a weak sense in The energy inequality ^ holds for any 0 < 
to <t. The uniform bound 

u{-, s)||i,oocis < 

holds with 

Koo = C J \uo{x)\'^dx 

and C a universal constant, independent of 6. If the initial vorticity ujq = 
V xuo is in \uJo\dx < 00 then it remains bounded in and moreover 



'R3 


u}{x,t)\dx— / \u{x,t)\^dx < / \uj{x,to)\dx— / \u{x,to)\^dx 


2 z/ 


'R3 


'R3 


2p 


'R3 


for all t >to. In addition the vorticity direction 

\uj{x,t)\ 

defined in the region {x||a;(a;, f)| > 0 } satisfies 


\uj{x,s)\\Vf{x,s)\'^ dxds < -u ^ 



0 ./{a:||a;(a;,s)|>0} 


'R3 


\uo{x)\'^dx. 


The proof of this result starts with the vorticity equation for the mollified 
equation: 


{dt [m] ■ V - A) a;* -f- eijkdj{[u]i){diUk) = 0. 


( 20 ) 




Here Cijk is the signature of the permutation (1, 2, 3) h->■ {i,j, k) and repeated 
indices are summed. Multiplying scalarly by ^ one obtains 

{dt+[u] ■ V-A) \uj\ + iy\uj\\V^\^ + Det{^,V[ui],diu) = 0 

where Det{a, b, c) is the determinant of the matrix formed by the three vec¬ 
tors a, b, c. Integrating in space and using the energy dissipation one can 
deduce the bounds for ut in and the bound on the direction For the 
bound on u in L°^ one uses the enstrophy differential inequality 

7 ^/ + u [ \Vuj{x,t)\^dx < 

jR3 JR3 




R3 


\uj{x, t)'\^dx 


'R3 


I Va;(a;, t)'^dx 


( 21 ) 


This is obtained from the vorticity equation (poD above by multiplication by 
cn, integration by parts, and use of bound (see (p7|) below) 


\ u \\ l ° 




( 22 ) 


Then one employs the idea of |p0 


one 


divides by (c^ + \u:{x^t)\^dxY^ 

and integrates in time using the energy principle, (c > 0 is a constant). One 
obtains a bound for 


/R3 


I Va;(a;, s) pda; > ds 


in terms of the initial data. The bound follows from interpolation (^21), 
the bound above, and the energy principle. We omit further details. 


4 Non-uniform bounds 

If M is a smooth solution of the Navier-Stokes equations and if 


'R3 


a;(x, t)\^dx ] dt < < oo 


(23) 
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then one can bound any derivative of u on (0, T] in terms of the initial data, 
viscosity, T and D (see [^, and references therein). Similarly, if one has a 
bound 


ll“(■.^)lli»(fc)* < B <co 


(24) 


then one can bound any derivative of u on (0, T] in terms of the initial data, 
viscosity, T and B. The quantities D and B have same dimensional count 
as viscosity (units of length squared per time). If one has a regularization 
that respects the energy dissipation and one has uniform bounds for the 
corresponding quantities then one can prove global existence of smooth so¬ 
lutions. If any of the two conditions is met then the solution is real analytic 
for positive times. Consider the mollihed equation (|^ at 5 > 0. Assum¬ 
ing for instance (|^) one obtains bounds for supj<'r \uj{x,t)\‘^dx and for 


Io^Irs \'^^{x,t)\‘^dxdt in terms of initial data, D and T, directly from (^Tj). 
The interpolation (^) then produces a bound for B. Vice-versa, if one has 
the assumption (|^) then one does not uses interpolation when one derives 
the enstrophy inequality from (^); rather, one integrates by parts to reveal 
u and one uses directly the assumption about ||m||l°o to deduce a uniform 
bound for the maximum enstrophy in the time interval [0,T]: 


sup / \uj{x,t)\‘^dx < S < oo (25) 

t<T Jr3 

This allows to bound D. In either case, the number S depends on the numbers 
D (respectively B) of assumptions (|^) (respectively (^^ ). By increasing £, 
if necessary, we may assume, without loss of generality the condition 

p = £^Tu-^ > 1. (26) 


This condition reflects the fact that we are not pursuing decay estimates. 
If no assumption is made then £ depends on h > 0. Once the enstrophy 
is bounded in time, higher derivatives are bounded using the Gevrey-class 


method of 28 


Theorem 3. Let 6 > 0. Consider solutions ofm) with initial data Uq E , 
cuo G Assume that one of the inequalities or (j^[) holds on the 

interval of time [0,T]. Then there exists a constant cq G (0,1) depending 
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only on the number p = p{£, u, T) of ^Bj), so that 


sup 

to<t<T 



e^^\^\\Q{f,t)\^df<2S 


holds for all 0 < to <T. Here A = ^/uT min cq}. If £ is uniform in 6 as 
5 —> 0, then the solution of the Navier-Stokes equations with initial data Uq 
is real-analytic and obeys the bound above. 


Using the fact that u = V x with divergence-free stream vector ip one 
sees easily that 

|2(0| < i|3(0| 

holds pointwise. It is elementary then to check that 










1 

4 


(27) 


holds for any positive A. Using the Fourier transform of ([TlD and the in¬ 
equality above one may follow closely the idea of |^. One considers the 
quantity 

y(t)^ [ 

JR3 


with V 
form 


^ and 0 < s < f < T and derives a differential inequality of the 


dy . (y . , - 

^<c(-) +c-y 


with absolute constant c. At f = s one has by construction y{s) < £. The 
differential inequality guarantees that y does not exceed 2£ on a time interval 
s <t < s + 2 cqT with Co a small non-dimensional constant (proportional to 
p~^ = u^T~^£~‘^). The time step 2coT is uniform because of the assumption 
that £ is hnite. One starts from s = 0. The differential inequality implies a 
non-trivial Gevrey-class bound for the second half of the hrst time interval 
[0,2coT]. One now sets s = cfT and makes another step of duration 2cqT. 
At each step the second half of the time interval yields a nontrivial bound, 
and because the initial point s is advanced by a half step one covers all of 
[cqT, T] . The pre-factor in the exponential bound is uniformly bounded below 


11 


by 2cqvT. If to < CqT we can advance 2to at a time, and obtain similarly a 
nniform bonnd on [to,T], with exponent pre-factor bonnded below by 2tQV. 
This resnlt implies that the velocity is real analytic. One may extend m to a 
complex domain z = x + iy hy setting 

u{x + iy,t) = {2 ti)~^ f 

iR3 

and, in view of (^7|) , the integral is absolntely convergent for \y\ < X, t > to. 
Note also that, at hxed 5 > 0 one has a hnite, (5-dependent bonnd on S, and 
therefore the solntions of the mollihed eqnations are real-analytic. Likewise, 
if one allows for 6 dependent bonnds, then the vortex method also can be 
shown to have real-analytic solntions. In this case, however, becanse the 
energy principle is not strong enongh one starts from the assnmption (^): 

Theorem 4. Let the initial vorticity loq belong to L^. Consider, ford > 0 the 
solution of ( fT^ ) and assume that holds on a time interval [0,T]. Then 
there exists a constant cq G (0,1) depending only on the number p = p(T, n, T) 
of so that 

snp [ <28 

to<t<T JR3 

holds for all 0 < to <T. Here A = y/i/T min cq}. If S is uniform in 6 as 
(5 —> 0, then the solution of the Navier-Stokes eguations ^ with initial data 
uio is real-analytic and obeys the bound above. 

The proof follows the same ideas as the proof of the corresponding resnlt 
for the mollihed eqnation ( 0 . 


5 Euler equations 

The three dimensional Enler eqnations 


dtu -\- u ■ Vu -I- Vp = 0, V ■ M = 0 


(28) 


are locally well-posed |^, |^. They conserve kinetic energy (if the solntions 


are smooth enongh [^, []^). Snch smooth solntions can be interpreted 
[jl| as geodesic paths on an inhnite dimensional gronp of transformations. 
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Despite of energy conservation, gradients of solntions may grow 
vorticity uj = V x u obeys 


m. The 


{dt + u ■ V) uj = uj ■ Vm. 


(29) 


Because of the quadratic nature of this equation and the fact that the strain 
matrix 


S = i{(V«)) + (V«)'} 


(30) 


is related to the vorticity by a linear classical singular Calderon-Zygmund 
integral, it was suggested that blow up of the vorticity might occur in 
hnite time. This problem is open, despite much research |5[|, [f^, HIT 


22 


29 


The blow up cannot occur unless the time integral of the maximum 


modulus of vorticity diverges 0. The vorticity magnitude obeys 

{dt + u ■ V) |a;| = a|a;|. 


(31) 


where the logarithmic material stretching rate a can be represented as 
a{x,t) = ^P.V. j D{y,i{x-y,t),i{x,t))\uj{x-y,t)\^. (32) 

Here 


D{y,i{,x-y,t),i{x,t)) = {y ■ i{x,t)) Det{y,i{x - y,t),i{x,t)) (33) 


and 




ui{x, t) 

uj(x,t)\' 


(34) 


In the two dimensional case ^ = (0, 0,1) and so a = 0. In three dimensions, 
if the vorticity direction is well-behaved locally in regions of high vortic¬ 
ity, then there is a geometric depletion of nonlinearity. More precisely, if 
^{x,t) X ^{x — y,t) vanishes in a quantitatively controlled fashion as ?/ —> 0, 
(for instance |^(a; — y,t) x ^{y,t)\ < k\y\), then a can be bounded in terms 
of less singular integrals (for instance in terms of velocity instead of vortic¬ 
ity). This observation suggests a correlation between vorticity growth and 
the geometry of vortex tubes. Such a correlation has been observed in nu¬ 
merical studies and was exploited to prove conditional results regarding blow 
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up for the Euler equations, 


, the Navier-Stokes equations [|^], and the 


quasi-geostrophic model [25|. The quasi-geostrophic model H, [17|, [23], [24 


40| is an example of an active scalar. Active scalars are advection-diffusion 
evolution equations for scalar quantities advected by an incompressible ve¬ 
locity they create: the velocity is obtained from the scalar using a hxed, 
time-independent formula: 


{dt + u- V — kA) q = 0 u = U[q]. 

The Euler equations themselves are an active vector system [IK 

(c?t -h M ■ V) A = 0, 

with 

u = 1E[A] 

the Weber formula [^ 

u{x,t) = P((VA)W). 


(35) 


(36) 


(37) 


(38) 


Here P is the Leray-Hodge projector on divergence-free functions, and 


v(x,t) = Uo(A(x,t)) 


(39) 


is a solution of 


[dt + u ■ V) V = Q (40) 

The initial data for A is the identity 

A(a:, 0)=a:, (41) 

and the initial datum for v is uq. Thus A{x,t) = a is the inverse of the La- 
grangian path a X {a, t)- The familiar Cauchy formula is in this language 


UJ 




1 



' dA 

’ dxi ’ dxj ) 


Here 


(42) 


C(x, i) =uio{A{x,t)) 


(43) 
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is the solution of 


(a* + u ■ V) C = 0 


(44) 


with initial datum ({x, 0) = uJo{x). 

One may use a near-identity approach to the incompressible Euler equa¬ 
tions: One solves (^) for short time, as long as VA — / is not too large. 
Then one stops and resets u = W[A\ in place of uq, sets A = x, adjusts the 
clock, and starts again. This approach allows one to interpret the condition 
PI of absence of blow up in terms of VA : if 

T 

\\VA{-,t)\\\oc,dt < oo 



then the solution of the Euler equations is smooth [0- 


6 Diffusive Lagrangian transformations 


The central object in the Lagrangian description of fluids is the Lagrangian 
path transformation a i—*• X{a,t)', x = X{a,t) represents the position at 
time t of the fluid particle that started at t = 0 from a. At time t = 0 
the transformation is the identity, X(a,0) = a. An Eulerian-Lagrangian 
formulation of the Navier-Stokes equations [|^ 

0 . 


parallels the active vector 


formulation of the Euler equations 


In order to unify the exposition we 


associate to a given divergence-free velocity u{x,t) the operator 

dt + u-V — uA = T^{u, V). 


(45) 


Using Ti,{u, V) we associate to any divergence-free, time dependent velocity 
held u a transformation x i—>■ A{x,t) that obeys 

rj,(n, V)A = 0, (46) 


with initial data 


A(x, 0) = X (47) 

Boundary conditions are imposed by considering the displacement vector 

£{x,t) = A{x,t) — X (48) 
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that joins the Eulerian position x to the the diffnsive label A. This is reqnired 
to vanish at inhnity, and one can think of A as being compnted by solving 

{dt + u- 'V — z/A) i + u = 0 (49) 

with initial data 

i{x,0) = 0. (50) 

If u is the solntion of the Enler eqnation and z/ = 0 in the eqnations (^, ^9]) 
then the map A is the inverse of the particle trajectory map a ^ x = X{a,t)- 
In the presence of viscosity this map obeys a diffnsive eqnation, departing 
thns from its conventional interpretation as inverse of particle trajectories. 
Nevertheless, continning the analogy with the inviscid sitnation, one nses 
the map x i—*• A{x, t) to pnll back the Lagrangian differentiation with respect 
to particle position and write it in Eulerian coordinates. This Eulerian- 
Lagrangian derivative is given by 

= Q*V (51) 

where 

Q{x,t) = {^A{x,t))~\ (52) 

that is 

Vf = Qjidj. (53) 

In the case z/ = 0 the invertibility of VA follows from incompressibility; in 
the diffusive case the determinant of VA does not remain identically equal 
to one as time passes. This imposes a constraint on the time of integration. 
We consider a small non-dimensional parameter g > 0 and work with the 
constraint 


sup sup |V£(a;,f)| < g. (54) 

0<i<r a;eR3 

With this constraint satisfied we can guarantee the invertibility of VA. In 
order to describe the dynamics and their relationship to the Eulerian dy¬ 
namics one needs to consider second derivatives of A. These influence the 
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dynamics because commutators between Eulerian-Lagrangian and Eulerian 
derivatives do not vanish, in general: 


[Vf.Vi] (55) 

The coefficients are given by 

= {Vf(d,£J}. (56) 


Note that 

= Q,,d,dkA^ = Vf{duA^) = 

These commutator coefficients are related to the Christoffel coefficients of 
the usual flat connection in computed using the change of variables a = 
A{x,t). With this change of variables, a straight line in x, x(s) = ms + b 
becomes the label path a{s) = A{x{s),t) and the geodesic equation ^ = 0 
becomes ^ = 0 with 

■pm _ /w 

ij — 

The simple geometry of is hidden behind a complicated transformation, 
but the transformation is the main object of study. The coefficients (but 
not u) enter the commutation relation between the Eulerian-Lagrangian label 
derivative and ri.(M, V): 

[T,(u,V)yf]^2i,CtA^i (57) 


This commutation relation is the viscous counterpart of the inviscid commu¬ 
tation of time and label derivatives. 

In the inviscid case the map A is the main active ingredient in the dynam¬ 
ics. The Weber formula (|3^) computes the velocity at time t directly 


from the gradient of A using a passively advected velocity v 
viscous case, the Weber formula 


In the 


u = P ((VA)*-?;) (58) 

can still be used but v{x,t) is no longer passive. Instead of (|i0| ) v obeys 

T^{u, V)n = 2i/C'Vn, (59) 
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that is 


r,(M,V)i;, = 2z/C'™afcT;^ (60) 

with initial data 


f (x, 0) = uo{x). 

The equations (|46D , (1^) together with the Weber formula 
to the Navier-Stokes equations [|T3: 


(61) 

are equivalent 


Theorem 5. Let A, v and u solve the system (^, ^). Then u obeys the 

incompressible Navier-Stokes equations, 


dtu + u ■ Vu — uAu + Vp = 0, V ■ M = 0. 


We describe now kinematic consequences of 


One starts with a 


smooth arbitrary incompressible velocity held u, one computes A using (^61), 
then one computes the inverse matrix Q = (Vv4)“^ (at least for a short 
time), and then Eulerian-Lagrangian derivatives and coefficients C. One 
then can evolve a vector v from an initial datum solving (^), and one can 
compute its Eulerian-Lagrangian curl C. = x v. The resulting equations 
r^ , |T5[] are summarized below. 


Theorem 6. Let u be an arbitrary divergence-free function and associate to 
it a map A solving and a vector field v solving ([5^). Then w defined by 

Wi = {diA^)vm (62) 

obeys the cotangent equation 

Vy{u,'W)w-\-{Wu)*w = t). (63) 

The Eulerian curl ofw,Vxw obeys the equation 

rj,(M, V) (V X tc) = (Vn)(V X w). (64) 

The Eulerian-Lagrangian curl ofv,( = x v, obeys 

V)C, = - 2i,Clja,Ci + (65) 
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The Eulerian curl of w, V xw and Eulerian-Lagrangian curl ofv,( = x v 
are related by the formula 


(Ve X w)g 



■ dA dA^\ 
’ dxi ’ dxj ) 


The determinant ofVA obeys 


r„(«, V) (log£>ei(VA)) = . 


( 66 ) 


(67) 


These considerations apply to arbitrary u without having to impose the 
equation of state (^8|). When z/ = 0 the cotangent equation (|63D appears in 
Hamiltonian formalisms [^], for the Euler equation in various gauges. 
The numerical merits of these have been analyzed critically H. When u is 
related to w by a hltered Weber formula, a gauge of the cotangent equation 
appears as a model of formally “averaged” Euler equation il and, in a 
viscous case, as a model of Reynolds’ equation P . 


Note that, when e = 0, the equation (|65|) is just the pure advection 
equation (5* + m • V) C = 0. If z/ > 0, C obeys a linear dissipative equation 
with coefficients Using Schwartz inequalities only, one obtains 


r..(«, V)|C(x,i)| +ii\VC{x,t)\ < 17i^|C(a:,«)| |C(a:,«)|' 


( 68 ) 


where 


\C{x,t)\^ 


C^.{x,t)C^,{x,t), \ax,t)\^ 


Cq{x,t)Cq{x,t) 


are squares of Euclidean norms. 
The evolution of the coefficients 


defined in (|56D is given by 


r„('t., V) (Gg) = -(5,7l’")vy4(%)) 
-(8,(«,))C,'” + 2 ^^,. a, (eg), 


(69) 


At time t = 0 the coefficients vanish, C]fi{x,0) = 0. Note that the linear 
equation (M) is identical to the nonlinear vorticity equation (^). It is the 
relation (|55D that decides whether or not we are solving the Navier-Stokes 
equation; if m is a solution of the Navier-Stokes equations then (pSD means 
u = Pic and consequently V x m = V x w implies that the Eulerian curl of 
w is the fluid’s vorticity. The formula (|66|) relating the vorticity to C, is then 
a viscous counterpart of the Cauchy formula. 
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Theorem 7. If u solves the Navier-Stokes equation, A solves (^^), and v 
solves ( \^ ) then the Eulerian curl of u, uj = V x u is related to the Eulerian- 
Lagrangian curl of v, ( = x v by the Cauchy formula 


00 


<? 


1 



■ dA 

’ dxi ’ dxj ) 


(70) 


Because of the linear algebra identity 


((VA)-^C), = {DetiyA)) 


-1 ^ 




Det 


C; 


dA dA 


dxi ’ dxj 


one has 


a; = {DetyA)) yA)-\. 
In two-dimensions become 

w=(Ce((V7l))C. 


(71) 


(72) 


reflecting the fact that, for z/ = 0, a; 
or (^) is the identity 


C, in that case. A consequence of (|70[) 


0 ^ ■ Vij = {DetyA)) (C ■ V^) 


(73) 


that generalizes the corresponding inviscid identity. Let us consider the ex¬ 
pression 


C{q,M) = {DetM)M-\ (74) 

dehned for any pair {q,M), where q G R^, M G GL(R^) are, respectively, 
a vector and an invertible matrix. This expression, underlying the Cauchy 
formula, is linear in q and quadratic in M, 

C{q,M)k = ^eijkDet{M,^i,M,j,q). (75) 

The quadratic expression in the right hand side is dehned for any matrix M. 
It is easy to check that 


C{q,MN)=C{C{q,M),N) 


(76) 
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and 


C{q,l) = q (77) 

hold, so C describes an action of GL(R^) in A third property follows 
from the explicit qnadratic expression ([75|) 

C(g, 1 + N) = {1 + Tr{N))q - Nq + C{q, N) (78) 


Here N is any matrix, and the meaning of C{q,N) is given by (^). If we 
consider, instead of vectors q and matrices M, vector valned fnnctions q{x, t) 
and matrix valned M{x, t) and nse the same formnla 

C{q, M){x, t) = C{q{x, t), M{x, t)) = {Det{M{x, t)) {M{x, t))~^ q{x, t) (79) 


then the properties ([76| , [77| , |78|) as well as 
Denote 

the total instantaneons energy dissipation rate, and 

K{t) = - [ \u{x,t)f dx 


obvionsly still hold. 


the total kinetic energy. We consider fnnctions u{x,t) that satisfy 


K{t) + / e{s)ds < Kq, 


( 80 ) 


and 


|^(*7 ^ 


( 81 ) 


For solntions of the Navier-Stokes eqnations and for solntions of the mollihed 
eqnation (^) the constants Kq, depend on the initial kinetic energy, vis¬ 
cosity and time only. For vortex methods, however, these constants depend 
on the cnt-off scale 5. The displacement £ satishes certain bonnds that follow 
from the bonnds above and (^91). 
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Theorem 8. Assume that the vector valued function i obeys ^ for t e 
[ti,T] and i{-,ti) = 0, for some ti > 0. Assume that the velocity u{x,t) is a 
divergence-free function that satisfies the bounds and G) on the time 
interval [0,T]. Then i satisfies the inequality 


||«(.,«)!U»< A-. 


(82) 


together with 


t)\^dx < 2KQ{t — tiY, 


(83) 


and 



\Vi{x, s)\‘^dxds < 


Koit-hf 


(84) 


ftl 



\'Vi{x,t)\‘^dxu / / \Ai{x, s)\'^dxds < Ci 


Ko{t-h) K^Eq 

V Z/2 


(85) 


Let us consider the analytic norms 




( 86 ) 


We will use p = 1 and p = 2. One can prove 

Theorem 9. Let the vector valued function i solve for t >ti > 0 with 
i{-,ti) = 0- Assume that a velocity u{x,t) defined for t > to , to < ti is a 
divergence-free function that satisfies 

sup ||M||{A,r,l} < Ur. (87) 

to<t<T 

Then there exists an absolute constant (a pure number) c such that 

||^(■,^)||{yl,r,l} < {t - ti)Ure''^ " (88) 


holds for t>ti. 


22 






(89) 


Note that (0) reads 

Note also that, if r < A then 


" ||Vo;p 


{A,A,2} 


II Vo; II {A,r,2} < 


-7T-r||<A^||{A,A,2} 

e(A — r) 


and 




(90) 


(91) 


hold. Combining Theorem 3 or Theorem 4 with the preceding result we 
obtain therefore 

Theorem 10. Consider solutions of the mollified equations o or of the 
vortex method ( \r^ ) associated to a filter Assume the initial data are 
divergence-free and belong to mq G V x mq = wq G Consider 
0 < f < T and let £ denote a bound for the enstrophy on the time interval 

10, TJ 

sup / \u}{x,t)'^dx < £. 

0<t<T Ar3 

Consider an arbitrary transience time 0 < to < T and the length scale 


A = -y/^min (to; Ciiy^£ 


with Cl a certain absolute constant. Then the velocity u obeys the bound & 
for any r < \, with Ur given by 


Ur = C2(A - r) 2^:2 


(92) 


and C 2 = 2a/ — . Consequently, for any ti > to, the solution i of the equation 


( UM) with initial data i{-,ti) = 0 obeys the bound and, for arbitrary 
ri < r < X, 


||V£(-,t)i|{A,ri,i} < e ^(r-ri) ^{t-ti)Ure' 

holds for t >ti. 


.(t-ti)ut 


(93) 
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These bounds depend on the cut-off scale S of the hlter only through the 
bound on enstrophy S. In particular, if the enstrophy is bounded uniformly 
for small 5 on a time interval then the results above apply for the Navier- 
Stokes equations on that time interval. If we measure length in units of a/ i/T 
and time in units of T then the enstrophy bound for u{x, t) = t) 

becomes 

sup ||Vu(-,s)||l2 < G 

0<s<l 

with G the non-dimensional number given by 

(94) 


Note that the number p of (|^) is just p = G^. In these units G is the only 
solution dependent parameter that we do not control. In terms of G and in 
these units, the dehnition of A becomes 


A 




= A = minjso; CiG }, 


(95) 


with So £ (0,1) arbitrary. The bound (|^ , with for 0 < r < A, r = a/ z/Tr, 
becomes 


Y — f/r = Uf < C 2 (A — r) 2 G. 

(96) 

For hxed sq and large G we may take A ~ G~^ so that if we 
(1 — 7 )A with 0 < 7 < 1/4 we deduce 

consider r = 

Uf < CsG^ 

(97) 

Choosing ri = (1 — j)r in (pSD we deduce from (PB|) 


l|V^(-,t)||{A;n,i} < 9 

(98) 

for to < ti < t < ti + tT with 


T = CigG~\ 

(99) 

Choosing r 2 = (1 — 7 )ri we deduce that 


\\VWi{;t)\\{A;r,p}<C,{uT)-'^G^g 

(100) 

holds on the same time interval. 
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7 Conclusions 


The viscous Navier-Stokes equations and their approximations can be de¬ 
scribed using diffusive, near-identity transformations. The velocity is ob¬ 
tained from the near-identity transformation using a Weber formula and a 
virtual velocity. The vorticity is obtained from the near-identity transforma¬ 
tion using a Cauchy formula and a virtual vorticity. The virtual velocity and 
the virtual vorticity obey diffusive equations, which reduce to passive advec- 
tion formally, if the viscosity is zero. Apart from being proportional to the 
viscosity, the coefficients of these diffusion equations involve second deriva¬ 
tives of the near identity transformation and are related to the Christoffel 
coefficients. If and when the near-identity transformation departs excessively 
from the identity, one resets the calculation. Lower bounds on the minimum 
time between two successive resettings are given in terms of the maximum 
enstrophy. 
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